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Abstract 

With a consistent definition of transverse-momentum dependent (TMD) light-cone wave 
function, we show that the amplitude for the process 7*7r° — > 7 can be factorized when the 
virtuality of the initial photon is large. In contrast to the collinear factorization in which the 
amplitude is factorized as a convolution of the standard light-cone wave function and a hard 
part, the TMD factorization yields a convolution of a TMD light-cone wave function, a soft 
factor and a hard part. We explicitly show that the TMD factorization holds at one loop level. 
It is expected that the factorization holds beyond one-loop level because the cancelation of soft 
divergences is on a diagram-by-diagram basis. We also show that the TMD factorization helps 
to resum large logarithms of type In x. 
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1. Introduction 

Factorization is an important tool to make predictions from QCD for an exclusive process if it 
contains both long- and short-distance effects, the latter involving one or more large momentum 
transfers, denoted generically as Q. It has been proposed long time ago that the amplitude of 
such a process can be factorized by an expansion of the amplitude in 1/Q, corresponding to an 
expansion of QCD operators in twist [TJ [2]. The leading term can be factorized as a convolution of 
a hard part and light-cone wave functions of hadrons. The light-cone wave functions are defined 
with QCD operators, and the hard part describes hard scattering of partons at short distances. 
In this factorization the transverse momenta of partons in parent hadrons are also expanded in 
the hard scattering part. At leading twist they are neglected. This results in that the light-cone 
wave functions depend only on the longitudinal momentum fractions carried by partons, not on 
transverse momenta of partons. This is the so-called collinear factorization. 

The collinear factorization has been also used for inclusive processes. At leading twist the 
transverse momenta of partons are also neglected, and the resulting standard parton distributions 
only depend on longitudinal momentum fractions. It is interesting to note that a transverse- 
momentum dependent (TMD) factorization can be made by introducing TMD parton distributions 
OHMS]- With a consistent definition, TMD parton distributions have some interesting relations to 
the standard Feynman distributions. For example, non-light-like gauge links shall be used to avoid 
some inconsistencies as discussed in the detail in [6]. The TMD factorization has many advantages. 
Beside those for different aspects of relevant physics, a prominent advantage is that it leads to the 
so-called CSS resummation formalism for resumming large logarithms appearing in the collinear 
factorization [31 [U [5]. This implies that the two factorization approaches are closely related to each 
other. 

In the collinear factorization for an exclusive process one usually has large contributions from 
end-point regions where the longitudinal momentum fraction x is small or close to 1, signified 
by large logarithms In 2 x in the perturbative expansion of the hard part. As such the expansion 
converges slowly and predictions can not be reliably made by taking few leading terms in the 
expansion only. It has been shown that those large logarithms can be resummed by using TMD- or 
kj- factorization [HJ. Similarly to inclusive cases, one introduces TMD light-cone wave functions 
(LCWF) by using non-light-like gauge links. Recently, such a factorization has attracted many 
attentions related to studies of exclusive I?-meson decays. For these decays, one can use collinear 
factorization to factorize short- and long-distance effects [Sj. Beside large logarithms, there are 
even end-point singularities in the factorization of certain decays. This makes the factorization at 
the end-point useless. In [TU] it has been suggested to use kr- or TMD factorization to resum large 
logarithms or to eliminate those end-point singularities. 

Comparing to TMD factorization for inclusive processes which has been studied extensively 
(see O HI [5] and references therein) , the TMD or /c^-factorization for exclusive processes have been 
not studied with consistent definitions of TMD LCWF's as much as that for inclusive processes. 
For a given process the amplitude in collinear factorization is factorized in terms of a hard part 
and LCWF's, latter representing all nonperturbative effects in the process. In [HJ [10] the TMD 
factorization can be taken as a " direct" generalization of the results of collinear factorization in the 
sense that the amplitude is factorized as a convolution only with a hard part and TMD LCWF's. 
From the study of TMD factorization for inclusive processes |3[ [U [5] one knows that the differential 
cross-section for a given process can be factorized as a convolution of a hard part, TMD parton 
distributions and a soft factor. The nonperturbative effects in the inclusive process are not only 
represented by the TMD parton distributions, but also by the soft factor. This is in contrast 
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to the collinear factorization, where the differential cross-section is a convolution of a hard part 
and standard parton distributions. Therefore, TMD factorization with just TMD LCWF's as 
nonperturbative objects in an exclusive process may be incomplete and some extra nonperturbative 
objects like soft factors may be needed to complete the factorization. It should be noted that the 
possible existence of soft factors in the factorization can change the scale evolutions which are used 
to resum large logarithms. One needs to examine the factorization in detail. 

In the study of TMD factorization for radiative leptonic decay of i?-meson|llj. it has been shown 
the decay amplitude takes the convolution form, with a hard part, the TMD LCWF of 5-meson, 
and a soft factor. The TMD LCWF of S-meson has been consistently defined and its relation to 
the standard LCWF has been studied in detail in [12]. From the result in [IT] the soft factor is 
really needed in order to make the factorization complete. A complete factorization means that 
the hard part does not contain any soft divergence like collinear- and infrared (I.R.) divergence. It 
has been shown that two factorizations for this simple case are equivalent and large logarithms can 
be conveniently resummed. 

In this work we study TMD factorization for the process 7*7r — > 7, where only one light hadron 
is involved. With the present study and that in [111 112] . we have examined TMD factorization 
for an exclusive process involving one hadron. This can be thought as a first step to the goal to 
examine TMD factorization for processes involving more than one light hadrons and to study the 
factorization for S-decays into light hadrons. A consistent definition of the TMD LCWF of it was 
first given in [7]. The TMD LCWF has been studied in[13|. where some interesting relations to 
the standard LCWF have been obtained. We will consider the case that the virtual photon has 
the negatively large virtuality —Q 2 (Q 2 > 0), which converts 7r into a real photon. Because of the 
large Q 2 , one can expand the amplitude in the inverse of Q 2 . We will show that the amplitude of 
7*7r — > 7 at the leading order of the inverse Q 2 can be factorized at one- loop level as the convolution 



In the above </>+ is the TMD LCWF, S is the soft factor and H is a coefficient function which can be 
safely calculated in perturbative QCD. The soft factor is defined with gauge links and it is needed 
to ensure H free from I.R. singularities. The correction to the factorized form is power-suppressed 
and is proportional to A 2 /Q 2 . A is any possible soft scale characterizing nonperturbative scale, like 
Aqcd, the mass of n etc.. This scale is about several hundreds MeV. Therefore, the correction to 
the factorized form is small with the large Q 2 . 

In the factorized form, the TMD light-cone wave function only depends on the specific hadron 
and is universal for all possible processes. The coefficient function H depends on the process but 
does not depend on the specific hadron. It represents a hard transition of a qq with the virtual 
photon into the real photon, in which the pairs move along one light-cone direction and has a 
small relative transverse momentum. The coefficient function H is obtained from the transition 
amplitude where the nonperturbative effects represented by soft divergences are subtracted. The 
soft-factor depends neither on the process nor on the specific hadron, it is completely determined 
by the infrared behavior of QCD. The TMD light-cone wave function and the soft factor are 
nonperturbative objects and can only be studied with nonperturbative methods of QCD, like Lattice 
QCD or QCD sum rule method. Phenomenologically, they can also be extracted from relevant 
processes. In the case of a pion, they can be extracted, e.g., from the process studied here, and 
those like it form factor, exclusive decays of S-mesons, etc... However, beside the process studied 
here, one needs first to show that the TMD factorization for other processes can be consistently 
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performed. In this work we will show that the proposed factorization for j*tt —* 7 holds at 
one-loop level and the cancelation of all soft divergences is on a diagram-by-diagram basis. The 
later is important to extend the factorization beyond one-loop level. We will also show that the 
factorization provides a simple way to resum large logarithms. 

Our paper is organized as the following: In Sect. 2 we give the definition of TMD LCWF and its 
one-loop results which will be used to show the factorization. In Sect. 3 we introduce our notation 
and show the TMD factorization at tree level. In Sect. 4 we show that the TMD factorization can 
be completed at one-loop level by introducing a soft factor. The soft factor is defined as products 
of gauge links. In Sect. 5 we show that possible large logarithms around end-point regions can be 
resummed with our method. Sect. 6 contains a summary. 



2. The Definition of TMD LCWF and Its One-loop Results 

We will use the light-cone coordinate system, in which a vector is expressed as a M = 
(a + , a~, aj_) = ((a + a 3 )/\/2, (a — a 3 )/\/2, a 1 , a 2 ) and a 2 = (a 1 ) 2 + (a 2 ) 2 . We introduce a vector 
= (u + , u~ , 0, 0) and the definition of TMD LCWF for ir° can be given as in the limit u + « u~: 

<p + (x,k ± ,(,v) = /^^ie i,!+ ^- i ^- £ M0|5(0)Lt(o t >,0)7 + 75i„(~,2)<!( Z )k (P))l z+ .(,, 

where q(x) is the light-quark field. L u is the gauge link in the direction u: 

L u (oo, z) = Pexp ^— ig s J dXu ■ G(Xu + z)^j . (3) 

The definition has been given in [7]. In the above, 7r° has the momentum P 11 = (P + , P~, 0, 0). 
The limit u + « u~ should be understood that we do not take the contributions proportional to 
any positive power of u + ju~ into account. The definition is gauge invariant in any non-singular 
gauge in which the gauge field is zero at space-time infinity. It has no light-cone singularities as 
we will show through our one-loop result, but it has an extra variable £ 2 . The evolution with the 
renormalization scale fi is simple: 

M 7^ = 27 9 + (a;,A;_L,C,At), (4) 

where 7^ is the anomalous dimension of the light quark field q in the axial gauge u ■ G = 0. 
In the above we do not specify if the light quark q is a u- or d-quark. The definition can be 
generalized to other light hadrons and there can be various TMD LCWF's. These TMD LCWF's 
of various hadrons have been classified in [TJ], and their asymptotic behavior with k± — > 00 has 
been derived [To]. 

To perform TMD factorization one needs to calculate the wave function with perturbative QCD, 
in which 7T° is replaced by a partonic state. We take the partonic state \q(k q ) , q(kq)) to replace ir 
in the definition, the momenta are given as 

k^ = (k+,k-,k q± ), k^ = (k+,kg,-k q± ), k+ = x P+, k+ = (l-x )P + = x P + . (5) 

These partons are on-shell. In this work we will keep the quark mass m for light quarks to regularize 
collinear singularities and give a small gluon mass A to regularize I.R. singularities. We will omit 
the /^-dependence in various quantities in this work and this dependence is always implied. 
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At tree-level, the wave function reads: 

0+ \x,k±,Q = S(x -x )S 2 (k ± - k q ±)<po, 4> = v(kg)-/ + -f 5 u(k q )/P + . (6) 

We will always write a quantity A as A = A^ +A^ -\ , where A^ and A^ stand for tree-level- 

and one-loop contribution respectively. At one-loop one can divide the corrections into a real part 
and a virtual part. The real part comes from contributions of Feynman diagrams given in Fig.l. 
The virtual part comes from contributions of Feynman diagrams given in Fig. 2., those contributions 
are proportional to the tree-level result. 




Figure 1: The real part of one-loop contribution to the TMD LCWF. The double lines represent 
the gauge links. 

The contributions from Fig. lb, Fig.lc and Fig. Id read: 

a s C F f / 1 \ x 1 1 + A 2 1 




2^2 j 
2 i \ 2 



(/)+(x,k ± ,C)\ lc = { 6{x - x Q ) ( — - — \ — ^ — ^— \s{x-xq) 5 1 ., ln g - 



x-x J + (l-x )A q 2 K u; gi + A 2 C 2 sg J™ 
<f>+(x,k±,()\id = ~ ~^f A 2 + q 2 5 ( x ~ g o)0o> (7) 
with 

= (1 - y) 2 m 2 + yA 2 + - y^x) 2 , y = — , q± = k ± - k q± . 

Xq 

A q = (1 - y) 2 m 2 + yA 2 + (k± - yk q± ) 2 , y = —, q± = k±-k q ±. (8) 

It should be noted that if we take u + = in the definition at the beginning there are light-cone 
singularities, represented by terms like 1/x and 1/(1 — x). With a nonzero, but small u + these 
singularities are regularized. All these contributions behave like kj 2 when k± becomes large. The 
contribution from Fig. la is too lengthy. But its calculation is standard because the gauge link is not 
involved. The contribution also behaves like fcj 2 when k± becomes large. The leading contribution 
for large k± is: 

2a 1/x 1 — x \ 

(f) + (x,k±,()\la = 37r 2 k 2 + A 2 \x~^ d ( X0 ~ ^ + 1 - Xo ( X ~ X °^J ^° 

+ "finite terms", (9) 

where we present only the U.V. divergent part explicitly and introduce a "cut-off" A, which also 
appear in the finite term so that the sum does not depend on it. All terms in the second line are 
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Figure 2: The virtual part of one-loop contribution to the TMD LCWF. The double lines represent 
the gauge links. 



proportional kj 4 when k± — * oo. Therefore, the one-loop correction of 4>+(x, k±, () is proportional 
to kj 2 when k± — ► oo, as expected from the general power counting rule |15j . 

The virtual part of the one-loop correction is from the Feynman diagrams in Fig. 2. Contribu- 
tions from each diagrams are: 



4> + (x, k±,()\2a 

4> + (x,k ± ,()\ 2b 

(f) + (x, k ± , C)| 2c 

(/>+(x,k±,()\2f 



<p+(x,k±,()\2d = <j)+\x,k±,C) 



6ir 



, u 2 , m 2 
m z X z 



<f>+(x, k ± ,()\ 2e = </>+(x,k±,C) ' ^ ln ^2> 



^\x,k ± ,() 



6tt 



21ni4 + 21n^+ln 2 ^-ln 2 ^ 



07T 



The complete one-loop contribution <j$/ is the sum of contributions from the 10 Feynman diagrams 
in Fig.l. and Fig. 2. 

From the above results one can derive some interesting relations between the TMD LCWF and 
standard LCWF and the evolution with £. The evolution with £ has been shown to be very useful 
for resumming large logarithms. Details about these can be found in [13|. 



2 In + 2 In 
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m 2 



X 2 



m 



X 2 



2~2 



+ br — T - In 

X z 



Cx_ 
X 2 



2tt z 

27T 2 

3 



+ 4 
+ 4 



(10) 



3. TMD Factorization for 7*7r° — > 7 at Tree-Level 

We consider the process: 

^°+7*^7 (11) 

where ir° carries the momentum P and the real photon momentum p. We take a coordinate system 
in which these momenta are: 

P" = (P+P- 0,0), p^ = (0,p- 0,0). (12) 

We will consider the case that the virtual photon has the large negative virtuality q 2 = (P — p) 2 = 
—2P + p~ = —Q 2 . The process can be described by matrix element {y(p, e*)|<77 M <?|7i"°(P)), which 
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can be parameterized with the form factor F(Q 2 ): 

( 7 (p,e*)|y/»g|7r°(P)) = ie^ e* v P pP(T F{Q 2 ). (13) 

For simplicity, we only consider one flavor quark q with electric charge 1. 

The form factor has been studied in experiments [16] . Theoretically it has been studied in the 
frame work of collinear factorization extensively [171 fl8| fl9j [20} [2~T| l22l 123] . In [22] the coefficient 
function in the collinear factorization has been calculated at one-loop, its two-loop result can be 
found in [23j . In [19] l20l [2l] an attempt has been made to take the transverse momentum into 
account in order to make predictions more reliable. In [21], following [2 [8] the resummation 
of large logarithms has been performed. However, in these studies, the issue whether the TMD 
factorization holds or not has not been investigated. The main goal of the paper is to address 
the issue, and we will not perform a phenomenological study with the factorization to make any 
numerical predictions. 

To find a factorized form for the form factor, we replace the hadronic state with the partonic 
state as before in calculation of the matrix element. The momenta of the partonic state is spec- 
ified in Eq.(5). In general one should neglect the transverse momentum k q ^ in the numerators 
in calculating F(Q 2 ), while one keeps the most important terms of k q ^ in the denominators [8]. 
Therefore, what neglected here is proportional to k 2 ± /Q 2 . However, at tree-level, the amplitude 
does not depend on k q ^. At tree-level there are two Feynman diagrams shown in Fig. 3. 




• ■ • 



Figure 3: Feynman diagrams of tree- level contributions to the partonic scattering. The black dot 
denotes the insertion of the electric current operator corresponding to the virtual photon. 



It is straightforward to obtain the tree-level result from Fig. 3: 



F(Q 2 



1 



|3a 



F(Q 2 



1 



1 



'Q 2 x ' * ^ A3b ru Q 2 (l-x ) ™QW 
Combining it with the tree-level result of the TMD LCWF, we can obtain a factorized form for the 

1 



(14) 



form factor: 



F{Q 2 



Q 2 



dxd k±<f)+(x, k±, C) 



1 1 
- + - 

X X 



(15) 



As discussed before, such a factorization may not hold beyond tree-level, and a possible soft 
factor should be implemented. We expect that the form factor in TMD factorization can be written 
as: 

1 



Q 2 
l 

Q 2 



dxdyd z k ± d 2 l ±( f> + (x, k±)S(y, l± t C)0(x + y)H(x + y) + 0(Q' 
®S®H + +0{Q~ A ), 



(16) 
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where we only give the dependence on convolution variables explicitly. The complete dependence of 
different objects on different variables will be given in the next section. With the above factorization 
one can easily find the tree- level results of the hard part H and the soft factor: 

H {0 \x) = - + -J- = 4, S^(y,l ± ) = 5(y)S 2 (l ± ). (17) 

X 1 — X XX 

It should be noted that at tree-level the hard part does not depend on transverse momenta. 
This is because we do not have an actual hadron in the final state of this simple case. When 
both initial- and final states contain hadrons, the hard part will depend on transverse momenta of 
partons entering hard scattering. 

4. TMD Factorization at One-Loop Level 

In this section we examine the TMD factorization at one-loop level and show that the soft 
factor is indeed needed. We will give its definition in terms of gauge links. At one-loop level, the 
factorization formula in the previous section can be written as: 

Q 2 jrW(Q2) = ® ® ff<°> + <^ 0) ® ® # (0) + </><£> 5<°> ® (18) 

If the factorization is right, will not contain any soft divergence and can be then calculated in 
perturbative theory. If all soft divergences in Q 2 F^(Q 2 ) are already contained in the first term, 
i.e., in <f>+ ® ® H^°\ then there is no need to introduce the soft factor. However, this is not 
the case. In this section we will study different factors introduced above. 

In deriving the factorization formula at tree-level, one should neglect the transverse momentum 
kqj_ i n the numerators in calculating F(Q 2 ), while one keeps the most important terms of k q ^ in 
the denominators [8j. As discussed in the last section, the hard part at tree level does not depend 
on the transverse momenta of partons. In calculating the one- loop contribution of F(Q 2 ) with the 
partonic state, we will therefore take the transverse momenta as small quantities and expand it 
in k 2 ± /Q 2 . We only take the leading order to extract the hard part. The higher orders are as 
power-suppressed contributions and they are neglected. 

4.1. The One-Loop Correction of F{Q 2 ) 

There are 12 Feynman diagrams for the one-loop correction of the form factor, 6 of them are 
given in Fig. 4. The other 6 diagrams are obtained from those in Fig. 4 by reversing the direction of 
the quark line. The diagrams in Fig. 4 represent the correction to Fig. 3a, and the other 6 diagrams 
for the correction to Fig. 3b. Two corrections are related each other by charge conjugation. Hence 
we will need to study how the contributions from Fig. 4 can be factorized. 



(a) 





('<) 



(0 



(</) 



(./) 



Figure 4: Feynman diagrams of the one-loop corrections to Fig. 3a. 
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The contributions except that of Fig.4a can be calculated in a straightforward way. They are: 



a R 



F(Q 2 )\ 4e = W)| 4/ =W)|3a~ 



_l n ii 7 _21n4-4 



F{Q 2 ) 



|4d 



F(Q 2 )U b 



F(Q 2 )Uc 



F(Q 2 )\ 3 a- 



F(Q 2 )\3a 



2a s 
3vr 

2a, 



In ^-r - In x + 1 

Q 2 



Q 2 l 

In + In — 7; + - In xq 



Q 2 



o7T Xq 



iln 2 x + 2Li 2 (-^) 

2 2 



1 

m — K in xo 

ttt/ 



+ 2 V 



In 



m 

Q 2 " 



2 + xq 
2x 



In xq 



(19) 



It is obvious that the contributions from Fig.4e and Fig.4f are already contained in the contribution 
of Fig.2a and Fig. 2d of the TMD LCWF, respectively. 

The contribution from Fig.4a is difficult to calculate and has a complicated expression. However, 
for the purpose of the factorization, we only need the contribution at the leading power of Q 2 . We 
have the result for the combination: 



3a 



-F(Q 2 )\sa 



2a s xo 
3tt 2 xq 



lnxo / d 2 k± 



1 



7T In Xq 



( Q 2 i N 

+7T In X I In -^j - 1 + - In Xq 



\ 2 + k\ 

(l + 0{Q' 2 )) , (20) 



In the above expression the I.R. divergent term In A 2 will be canceled from the integration of k±. 
The expression contain an U.V. divergence. This divergence comes from TMD LCWF of Fig. la 
and needs to be subtracted. 

Since the contributions from Fig.4e and Fig.4f are already contained in the corresponding contri- 
butions of TMD LCWF, i.e., they are already subtracted for the hard part, we have only collinear 
divergences in F(Q 2 )\^i ) and F(Q 2 )\4 C , signalled by mm, as soft divergences which need to be 
subtracted. 



4.2. The Convolution ^ ® ® 

In this section we present the results for the convolution <p^) <g) £K°) <g) . We will see that 
the convolution contains the same collinear divergences as those in F^ l \Q 2 ). Therefore, with our 
TMD LCWF the collinear singularities are correctly factorized. Moreover, the collinear divergence 
in a given diagram for F^ l \Q 2 ) appears exactly in the same class of diagrams for TMD LCWF. 
This fact can make it possible to extend the factorization of collinear divergences beyond one-loop 
level. We will also see that there are still some I.R. divergences in the convolution. Therefore, a 
soft factor is needed to subtract them. 

We find that it is convenient to present the results in the following combinations: 



lib + ' 



._,..„ .-^+|2c]®5 (0) ®^ (0) ^ 



1 f d 2 q ± 1 



vr ql+X 2 q 2 ± +X 2 



(■2 2 



9 




2Li(-^ 




(21) 



Comparing the above with the one-loop correction to the form factor in Eq. (|19p . we see that the 
collinear singularity in F(<3 2 )|4b and F(Q 2 )|4 C is exactly reproduced in Wb and W c , respectively. 
But, in the convolution <fi+ (g) ® there are also I.R. singularities shown as In A 2 and 



I.R. integrals in Eq. (|2ip . It should be noted that these integrals are also U.V. divergent. These 
U.V. divergences are closely combined with certain I.R. divergences. Once these I.R. divergences 
are subtracted, one can expect that the U.V. divergences are subtracted too. Because of these 
divergences, one needs to introduce the soft factor so that the factorization can be completed, i.e., 
the hard part H^ 1 ' subtracted from Eq. (|18p is free from any soft divergence. 

4.3. The Soft Factor 

As discussed in the above, the I.R. divergences which need to be subtracted are from TMD 
LCWF. In general these divergences can be analyzed by applying an eikonal approximation. Using 
the approximation for those diagrams in Fig.l. one can find that the I.R. divergences in these 
diagrams can be reproduced by replacing the quark lines with eikonal lines. The corresponding 
diagrams just represent the contributions to the product of the gauge links: 



1. 



Tr(0|T \lU0, -00)4(00,0)^(00, z)L vl (z, -oo) 



3 

Pexp 



Si(z ,z±,vi,v 2 ) 

L Vl 2 (z, -oo) = Pexp^-ig s J d\v\^ ■ G(Xv^2 + z) j , 
S i (y,q±,ClX2) = P + J dz-d 2 z ± e iyP+z -- t ^S 4 (z-,z ± ), 
the two vectors v\ j2 and the parameters £i 2 are defined as : 

2v7 2 (P + ) 2 



-1,2 



K2.»i,2.".0). v t2» v ia> Ci,: 



v l,2 



(22) 



(23) 



these vectors represent the trajectories of q and q in Fig.l, respectively. The limit vf 2 » v i2 
should be understood that we neglect all terms with a positive power of £1,2- The ordering of the 
limits is relevant for calculating S4. We will take the ordering that the limit » is taken 
before taking the limit v 2 » v 2 ■ 

We have calculated S4 in perturbation theory. At tree-level we have 



sf\y,q±XiX2) = 5(y)5' 2 (q ± ), S^'(z ,z±,v 1 ,v 2 ) 



(0), 



1. 



(24) 



At one-loop level, the contributions to £4 can be represented by those diagrams in Fig.l and Fig. 2, 
where one replaces the quark line with the gauge link L Vl and the antiquark line with the gauge 
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link L V2 . The contributions from Fig.l are: 



'S'4(y,^L,Cl,C2)|la 

S4(y,q±,Ci,(2)\ib 
5 , 4(y,g±,Ci,C2)|ic 
5 , 4(y,g±,Ci,C2)|id 



2a s 
3tt 2 

2a s 
3vr 2 

2a s 
3vr 2 



0(-y) 



+ 



</(<?! + A 2 + dV) y(9i + A 2 + C 2 V) 

0(-y) /1\ 1„, , 1 



In. 



2a 5 6(y) 
3vr 2 g 2 +A 2 " 



(25) 



The ordering of the limits mentioned above only affects the contribution S^\x a . If we change 
the ordering, S4I10 changes its sign. The above expressions, when they are used to calculate the 
convolution ® fiK 1 ) (8) , will be taken as distributions for — xq < y < 1 — xq. This can be 
seen from the convolution 



0f ® 5« ® i? (0) = </>o / dj/dgi-SW (y, g ± )0(* o + y)H® (x + y). 



(26) 



In Eq. ([25|) we have already take the limit ( — > 00, i.e., u + << u~ by taking them as distributions 
for — x < y < 1 — aro- 

The contributions from the diagrams corresponding to those in Fig. 2 are: 



Si(y,q±, C1X2) 



\2a 



Si(y,q±,(i,(2)\2b = Si(y,q±, (i,(2)\ 



2,1 



5 , 4(y,g±,Ci,C2)|2e = S\ '(y,q±, (1,(2) ' ^ ln ^> 



'S'4(y,^L,Cl,C2)|2 C 

S4(y,q±,(i,(2)\2f 



a. 



II 



-5'| 0) (y , q± , Ci , C2) ^ In ^ ln ^2 > 



c 2 

2 a2 



-£4° 9±> ft, (2)7^ In ^ In ^ 



(27) 



S4 has the charge conjugation symmetry £4(2/, <fj_, (1, (2) = 5*4(— y, — q±, (2, Ci)- We can use the 
symmetry to simplify our definition of the soft factor by taking the limit (2 —> (1 = (v 

If we identify the soft factor S with S4 as S(z~ , z±) = [S&(z~, z±, v\, V2)]~ l , the I.R. divergences 
in (f>+(x, k±, () are indeed subtracted if we integrate over x and k± and extend the integration region 
of x from —00 to 00. However, in these convolutions (f>+ <g> <8> and ® <8> H^°\ 
the integrands are not or (ff^ alone, but multiplied with and the integration region of 

x is restricted from to 1. This results in that the I.R. singularities in <fi+* ® are not 



exactly reproduced in cp\ ® (g> H^°\ and there are end-point sing ularities in $ ' <g> <g> 
when y approaches to — xq or 1 — xq. The problem is related that the proposed soft factor does not 
vanish at the end points, while the TMD LCWF does. Therefore the convolution with S(z~ , z±) = 
[S^z" , z±, v\, V2)}~ 1 is ill defined and the proposed soft factor is not right for the case. 

To overcome the problem, one can construct quantities from S4 and their contributions have a 
one-to-one correspondence to the one- loop contributions to (f>+(x, k±, () in Eq.([8]) and Eq.([9]). The 
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correspondence means that the I.R. divergences in Eq. ([8]) and the U.V. divergences in Eq.© are 
exactly reproduced. For this we first define: 



l-e(-y)- 



y 



y-x 



y 



y + l — x 



1 



lim [S±(q + ,q±, (1,(2) + S^(q + ,q±, (2,(1)] ■ 

£ C2^Cl=Cu 



(28) 



The contributions from Fig.l to S\ can be easily obtained from those of S4. We have: 
•Si(y,x,q±,Cv)\ia = 0, 



Si(y,x,q±,( v )\ib 



2a, 



3^ 2 P+ 



9(x - x) 



qj_ + A 2 + Q(x - x ) 2 x \x-x J + 



1 



In 



2^2 




C 2 x 



-\ 5{X - XQ) ^\x 2 ^ql + X 2 



Si(y,x,q±,Cv) 



2a, 



lc 



3vr 2 i^ 



1 



^ + A 2 + ( 2 {x - x J" x \x-x / + 

'2^2 



if* N 1 , C ^0 

+ 2^- Xo) ^T^ ln ^iTA^ 



Si{y,x,q±,(v)\id 



2a s 5(x — xq) 
3tt 2 P+ g 2 +A 2 



(29) 



where we have taken y = x — xq because of Eq. (|26p . The contributions from Fig. 2 to S\ are the same 
as those of S4 by taking £2 = Ci = Cv- Comparing Eq. ([29|) with Eq.©, we can see that the I.R. 
singularities in <f>+\i c and <p+\id are exactly reproduced <Si|i c and respectively. 

Next we define 



S2(y,x,q ± ,(v) 



y 

1 



l-O(-y)- 



y 



y-x 



0{y) 



y 



y + 1 — x 



- lim [S4:(q + , q±, (1, (2) ~ S<i(q + ,q±, (2,(1)] ■ 
<Z Ca— »Ci=C» 



(30) 



One easily finds that 1S2 only receives a nonzero contribution from Fig. la: 
S 2 



S2(y,x,q±,C 



v la 



<?2 1 la, 

2a s 
3vr 2 P+ 



8(x - x) 



+ 



9{X - Xq) 



1-x 



q\ + A 2 + Q{x ~ x o) 2 x o q\ + A 2 + Q(x - Xq) 2 1-Xq 



,(31) 



where we have taken y = x — xq again. It is noted that the U.V. divergence in </>+|i a in Eq.© is 
reproduced by <S2|ia- 

There are extra I.R. divergences introduced by the contributions from Fig. 2a and Fig. 2d for S4. 
These can be subtracted by introducing 



S 2 (Ci,C 2 ) = ^Tr(0|T[4 2 (0,-oo)L Ol (0, -00)] |0). 



(32) 
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Figure 5: One- loop contribution to S^- The double lines represent the two gauge links. One is for 
L Vl , the other one is for . 

As for S4, the same limits about the vectors v\ and i?2 are taken. At leading order = 1. At 
one-loop level, the contributions are from the diagrams given in Fig. 5. They are: 

$2 1 5a 



2\5b 



a s , ^ 2 , Ci 2 

3^ ln ^ ln cr 

CI _ ° s 1 ^ 
^2 5c — -5— in "To ' 



(33) 



In S\ we have taken the limit £2 — ► Ci = Cv Before the limit there is an I.R. singularity in 5i a . 
This singularity is exactly the same as S^lsa- We take the limit C2 — ► Ci = Cv for 52 and define: 



S 2s (Cv) = A lim i(S 2 (Ci,C2)+5 2 (C2,Ci))- 



(34) 



The I.R. singularities in Fig. 2a and Fig. 2d for S4, hence for S\ are reproduced by S^s- It should 
be noted that the limit (2 — ► Ci f° r is nontrivial. If we take £1 = £2 at the beginning, we will 
always have S2 = 1- 

Finally we can define the correct soft factor as: 



5 2 



5(z ,^±,x,Ct,) = 

%,ar,9l,C«) = P + J dz-d 2 z ± e^ z --^^S(z-,z ± ,x,(v). 
With this we will show that the extracted hard part is free from any soft divergence. 

4.4. The Convolution ^ <g> S"W ® # (0) and The Result for ffW 

With the soft factor we define the following combinations of convolution, 



2a. 
3ttxo 

2a. 

37TXQ 



(35) 



— + - In 2 
12 4 



1 A 2 C 2 1 . A 2 1 /■ A ± 
— In — ^ in -tt H — m — r H — 

2 /x 2 C« 2 n 2 2 

(0) 



gj + A 2 

vr + A 2 C 2 ^o 



In 



vr 2 1 1 2 A 2 
h - hr — r-^ 

12 4 C 2 ^ 



A 2 C 2 



2 1 V ln c 2 + 5 1 V " M 2Li( "? + ln Xo ln CM 



A 2 
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In 



XQ 



d 2 



1 



XQ 



7T 




H. 



?1 + a 2 

(0) _ 2a 



3<i 



In xq 



f d 2 q± 1 

/ 7T (? 2 +A 2 

dV 1 

37TX J tt qj_+ \ 2 



111 



3a 
(i 2 



q± 



XQ 



vr g 2 + A 2 



Xn 1 
+ 2 In — - — 

x x 



C 2 4 



2L* 2 (-^)-lnx ln%° 



(36) 



The convolution 



,(0) 



5(1) #(o) ig the sum Ua+ u b + u c + u d Comparing Eq.flU]) and Eq.flM 



one can see that the I.R. divergence and U.V. divergence in W&, W c and Wd is canceled by W&, W c 
and Wrf, respectively. The sum X)j=6cd + W»] contains only collinear singularities which are the 
same m FW(Q 2 ). The U.V. diver gence in Eq. (l20p is also canceled by that in U a - Hence, with the 
soft factor the hard part H will be infra-red finite. It is clear that the soft divergences are canceled 
or subtracted in a diagram-by-diagram basis, e.g., the collinear singularity in Fig. 4b for the form 
factor is subtracted by that in Fig. lb and Fig.lc for the TMD LCWF, and the I.R. singularity 
in in Fig. lb and Fig.lc for the TMD LCWF is canceled by that in in Fig. lb and Fig.lc for the 
soft factor. The contribution to the TMD LCWF from Fig. 2b is completely subtracted by the 
contribution from Fig. 2b to the soft factor. Therefore, after the subtraction and cancelation the 
contribution from Fig. 4b to the hard part is free from any soft divergence. 
We now can subtract the contribution from Fig. 4 to get H^: 



■le 



H(x 
H{x)\ 4d 

H{x)\ Ab 

H(x)\ Ac 



= H(x)U f = 0, 

,2 



-a s 
3irx 

3irx 
a 



hi — ^ — In £ + 1 



f 7T 2 O 2 1 C 2 

|-2 + y+ ln^-ln, + -lnS f 



c 2 c 2 

In -^r + In + 4 In x 

fj, 2 fj, 2 




-ling 

2 C 2 



ln^ +ln^| +41nx 



H{x)U a = 



2a s 
3irx 



, In a; 



Clx 2 



2 (l„, x x - I 

+ — m „ in x + 3— in x > , 
x Q z x 



1 



ln ^ + 2 -2 lnI 



- 2Li 2 (-~) + 2x hi- 
re x 



(37) 



Combining the results from diagrams obtained by reversing the direction of quark lines in Fig. 4 
' *° obtain the complete contribution to H^: 



wc 



2a, 1 



XX 



1, Q 2 , <5 2 , i , C 2 A C 2 . C V 2 V , C 2 w 



c _ c _ 

-2 In —77 ln(xx) + In — ^ (x\nx + xln x) 
Q 2 Q 2 



+ 



3irxx 



—5 — x In 2 x 



-x In 2 x + xlnx + xlnx — 2 ln(xx) + 4xLi2(x) + 4xLi2(x) 



(38) 
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The factorization formula with complete parameter-dependence now reads: 

F (Q 2 ) = Q2 J dxdyd 2 k_ L d 2 l_ L (f) + (x,k_ L ,C,fi)S(y,x,l±,C,Cv,fJ')0(x + y) 

■H(x + y,Q 2 ,C,(v,v) + 0(Q- i ), (39) 

with 

H(x, Q 2 ,C, C„ /x) = H^{x, Q 2 ,(, ( v ,») + H^(x, Q 2 , (, ( v , = 4 + H^(x, Q 2 , (, (40) 

xx 

from Eq. (|38p it is clear that the hard part H does not contain any soft divergence and we complete 
the TMD factorization at one-loop level. Since the soft divergences are canceled in a diagram-by 
diagram basis, the factorization can be argued to hold beyond one-loop level. Prom our results in 
this section, it is clear that a soft factor is needed so that the TMD factorization with the TMD 
LCWF defined in Eq.([2]) can be completed. 



5. Resummation of Large Logarithms 

We note that our factorization formula can be simplified by introducing <j>+(x, k±, £, ( v , fj,): 

dz- d 2 z ± ik+z 



v z —izj_-k 



With 



2vr (2vr) 2 

•(0| 9 -(0)4(oo,0) 7 + 75^(oo,z) ( 7(^K (P))| 2+= o. 
the factorization formula reads: 



F(Q 2 



1 

Q 2 



dxd 2 k ± 4> + (x,k ± ,(,(v,fJ')H(x,Q 2 ,(,(v,n) + 0(Q 



-4\ 



(41) 



(42) 



In principle one can re-define the TMD LCWF as that given in Eq. (|41|) . With the redefined TMD 
LCWF the factorization takes a simple form. However, we do not intend to do so because it is 
unclear if the same re-definition is useful for other processes whose TMD factorization have been 
not studied in detail yet. We will call <f>+ as a modified TMD LCWF in this work. It should be 
kept in mind that the soft factor is there. 

With the one-loop result for every factor, we can work out the evolution equations for the hard 
part, 



M ^lntf(x,Q 2 ,C 2 ,C,V) 
C-^lnH(x,Q 2 ,( 2 ,(lv 2 ) 
C v J^lnH(x,Q 2 ,( 2 ,(lv 2 ) 



2a s 
' 3vr 



2a s 
3vr 

2a s 
3vr 



t 2 l 

l + 21n^ +0{a 2 s ) 
C 2 

-2 + 21n-^ + 21n(xx) 



+ 0(a 2 ) 



-2 In 



V 2 



4 ln(xx) + 2 [x In x + x In x) 



+ 0(a 2 ). (43) 



The corresponding evolutions of 4> + reads: 
d 

W+(^ fc -L>C 2 >C 2 ,^ 2 ) : 

d 

C^in<Mx,fc_L,c 2 ,C 2 ,^ 2 ) : 



2as 
3vr 

2a s 
3vr 



l + 21n 



C 2 



+ 0(a 2 ) 



C 2 

2 -21n-^ -2ln(xx) 
H 1 



+ 0(a 2 ), 



(44) 
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the evolution with does not take a simple form, it takes a convolution form: 



Cv -Sr 4>+ (% ,k± , ( 2 , C p 2 ) 



with 



C v {x,y) 



2a s 
3?r 



5(x-y) 2 In 



CvVV 
P 2 



dyC v (x,y,(v,p)4>+(y,k±,( ,C„,M ) 



i 0(y ~ x)- - 6(x - y)^ 

x-y \ y y 



(45) 



+2[9{y- x)- + 9{x-y)- 

y y 



(46) 



where the last line comes from the box diagram. The +-prescription is for x. It should be noted 
that: 



dx^ 

o xx 



1 



x-y 

dx^ 

o xx 



6{y-x 



x 



y 



(x - y)'- 



V 



6(y — x) — h 6(x — y) — 

y y 



1 

yy 
l 

yy 



[2 + ln(yy)] . 
(y In y + y In y) . 



(47) 



With the above equations one can easily check that the form factor does not depend on (i, £ and 
( v , as expected. 

In Eq. (l38j) there are different possibly large logarithms. Beside those resulted from different 
scales like £ 2 , (%■■■ e ^ c -) there are large logarithms around the end points of x, like ln 2 x and lnx if 
x goes to 0, and ln 2 (l — x) and ln(l — a;) if x goes to 1. These logarithms also appear in the collinear 
factorization. They make the perturbative expansion of H unreliable and must be resummed. With 
the TMD factorization, the resummations can be accomplished in a simple way. 

The result for H can be simplified by noting that the modified TMD LCWF is symmetric under 
the exchange x <-> 1 — x. It is convenient to introduce a parameter p to relate £ 2 and C, 2 : 



2„,2 



c 2 



U V 



4(v ■ u) 2 

The factorization in Eq. (|42p now reads: 

1 



= C/p 2 , 



1 v~<~u 



V u n 



Q 2 



dxd k ± (f> + (x, k ± ,(, P, p)H(x, Q, (, p, p) + 0(Q ). 



(48) 



(49) 



and 



2a s 
3irxx 



1 Q 2 Q 2 1 o 2 2 C 2 
_ln_ + l n ___l n p + ln p l n _ 



/i- 



c 



+ 



2 In — (In x — x In x) + 2 In p (2 In x — x In x) 

x In 2 x + xhix — 21nx + AxL\2 (x) 

3irxx 2 



(50) 
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where we have used the property (f>+(x, k±, C, p, p) = </>+(l — x, k±, (, p, p) to simplify the hard part 
H . The evolution of H with various scales reads: 



d 

P-q^^H(x,Q,C,P,p) 



2a, 
' 3vr 



l + 21np 2 +0(a 



C,-^- \uH(x, Q, (, p, p) = 1 — In/? 2 + 21nx — 2xlnx 



p-7^ln-H(x,Q,(,P,p) = ^ 



C 

In p 2 + In — + 4 In x — 2x In x 

P 2 



+ 0(a 2 s ) 

+ 0(a 2 s ) 



(51) 



To avoid large logarithms in the wave function, we first choose appropriate scales to eliminate 
them ^ 

F ^ 2 ) = 7^2 i im n / dx4> + (x,( L ,PL,PL)H(x,Q,(L,PL,PL) + 0(Q- 4 ). (52) 

If we use the result in Eq. ([50j) . there will be large logarithms in H(x,Q,(l, Pl, Pl), especially 
terms like In 2 x, lnx when x — > 0. It is possible to resum those large logarithms terms by using the 
evolutions in Eq. (j51~j) to express H(x,Q,(l, Pl, Pl) as: 



H(x, Q, Cl,Pl,Pl) = e s H(x, Q, (h,Ph,Ph)- 



(53) 



This represents the resummed form for the hard part. By taking appreciated scales Qh, Ph and ph 
so that H(x,QXh, Ph, Ph) does not contain any large logarithms. It should be noted that lnp^ 
is not a large logarithms. We first use the /x-evolution to evolute the hard part to a high scale pn 
and take pn = Q. As the next we use the {^-evolution to evolute H to a high scale Qh = Q- After 
these the hard part reads: 



H{x,Q,(l,Pl,Pl) 



H(x,Q,(h = Q,Pl,Q) 





r dp 


exp < 






I Kh P 




{ 2a s (Q) r 


• exp 


\ 3vr i 




' 2a s (Q) 


XX 


3vr 



(l + 2 In pi 



1 — In p L + 2x In x 
1 



ln^ 

CI 



H(x,Q,( h ,Pl,Q), 



In 2 p\ + 2 In p L (2 In x — x In x) 



x In 2 x + xlnx — 21nx + 4xLi2 (x) 



(54) 



Now the hard part only contains possibly large logarithms In x and lnx when x goes to zero. To 
resum these logarithms we note that 



loo o o 1 

- — In p + 4 In p In x — In x — 2 In x = — 
2 H H 2 



In 2 C - In 2 P \ 



1 

V 



p 2 H = p 2 H (x)=x^e-V^. 
Using the p-evolution we have: 
H{x,Q,(l,Pl,Pl) = exp(-S) • H(x,Q,Q,p H ,Q), 



(55) 



2a s (Q) 
3vr 



c 2 

Cl 



1 — In p L + 2x In x In + ln z x + 2 In a; — - — 2x Inx In p H 
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H(x,Q,Q,p H ,Q) = 4U + 



+— In 2 p\ — lnp|(41nx — 2x lnx)l — / 
2 J J fi 

2a s (Q) 



*3 <i// 2a s (//) 

/i 37T 



(l + 21np|) 



3tt 



— — 2x In p\ In x + x In 2 x 



+x In x + 4xLi2 (x) — — 



(56) 



Now there is no large logarithms in H(x,Q,Q, pn,Q). It is instructive to compare the result of 
H(x,Q,£l, PL, P>l) with and without the resummation for x — > 0. Without the resummation the 
result for x — > reads: 



m n > \ -i 2u s {pl) . 2 
xH(x,Q,C L ,p L ,p L ) ~ 1 In x, 

07T 



(57) 



therefore the one- loop correction can be very large and is divergent. After the resummation the 
result for x — > reads: 

2a s (Q) , 

'" " (58) 



xH(x,Q,( L ,p L ,p L ) ~ x " S*" lna: , 



which will approaches to faster than any positive power of x. 



6. Summary 

In this work we start with a consistent definition of TMD LCWF's to explicitly show that the 
TMD factorization for the process 7*7r° — > 7 can be made at one- loop level, where a soft factor 
must be implemented. The soft factor is defined as a product of gauge links. It makes the hard 
part free from I.R. singularities. In contrast to the corresponding collinear factorization, where 
the amplitude is factorized as a convolution with the standard light-cone wave function and a 
hard part, the amplitude in the TMD factorization is factorized as a convolution of a TMD light- 
cone wave function, a soft factor and a hard part. From our results, the cancelation of all soft 
divergences happens in a diagram-by-diagram way, hence it is possible to extend the factorization 
beyond one-loop level. 

With the factorization, the hard part can be safely calculated in a perturbative expansion in 
a s . But in general results from perturbative expansions have large logarithms, which can make the 
results unreliable. A resummation of these large logarithms is needed. With the TMD factorization 
we show that the resummation in our case can be done simply. It should be mentioned that it is 
possible to show that the two factorization are equivalent, as the study of TMD factorization for 
the radiative leptonic decay of 5-meson shows [TT]. 

After the detailed examination of TMD factorization for the simple cases involving one hadron 
only, we can proceed to study more complicated cases involving more hadrons, particularly 5-decays 
into light hadrons. 
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